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We set up a forward-backward path integral for a point particle in a bath of photons to derive a 
master equation for the density matrix which describes electromagnetic dissipation and decoherence. 
As an application, we recalculate the Weisskopf-Wigner formula for the natural line width of an 
atomic state at zero temperature and find, in addition, the temperature broadening caused by the 
decoherence term. 



I. INTRODUCTION 



The time evolution of a quantum-mechanical density matrix p(x +Q , x_ a ; t a ) of a particle coupled to an external 
electromagnetic vector potential A(x, t) is determined by a forward-backward path integral [[l] 



(x +6 ,<b|x +a ,i a )(x_ 6 ,i b |x_ a ,< a )* = f7(x +b ,x_ 6 ,t 6 |x +a ,x_ Q ,i a ) 

ft b r 



2?x + 2?x_ exp < — , 
K Jt 



— (4 - x 2 _) - y(x+) + V(x_) - ; x + A(x + , t) + ^x_A(x_,<) 



(1.1) 



where x + (t) and x_(t) are two fluctuating paths connecting the initial and final points x +a and x +b , and x_ a and 
x_b, respectively. In terms of this expression, the density matrix p(x +b , x_ b ; i b ) at a time i b is found from that at an 
earlier time t a by the integral 



p(x +h , x_ b ; t b ) = J dx +a dx„ a U(x +b , x_ fc , i b |x +a , x„ a , t a )p(x +a , x_ a ; t a ). 



The vector potential A(x, t) is a superposition of oscillators Xk(i) of frequency f^k = c|k| in a volume V: 

s 4kx f d 3 kV 



A(x,t) = 5^Ck(x)X k (t), c k 



(2tt) 3 



(1.2) 



(1.3) 



These oscillators are assumed to be in equilibrium at a finite temperature T, where we shall write their time-ordered 
correlation functions as G£ k ,(i,f) = (T li k , (*')) = $wG Qk {t,t') = ~ WW /k 2 )Gn k {t, t'), the transverse 

Kronecker symbol resulting from the sum over the two polarization vectors Ylh=± e *(^-' h)e J *()s., 'O °f the vector 
potential A(x, f). For a single oscillator of frequency f2, one has for t > t': 



1 h cosh ^ -»(*-*')] 

G ,(M') = 5 [*,(*, + c„(t,0] = ^ 2 — m , t > t> 

sinh — - — 



(1.4) 



which is the analytic continuation of the periodic imaginary-time Green function to r = it. The decompostion into 
Aa(t,t') and Ca(t,t') distinguishes real and imaginary parts, which are commutator and anticommutator functions 
of the oscillator at temperature T: Cn(t,t') = ([X(t), X(t')]) T and A n (t,t') = ([X(t), X(t')]) T , respectively. The 
thermal average of the probability (1.1) is then given by the forward-backward path integral 



?7(x +b ,x_ b ,i b |x +a ,x_ a ,i a ) = J Vx + {t) J Px_(t) 



x exp < — / dt 



M 



(x* -x 2 _)-(F(x + )-V(x_)) 



A FV \x + , X J 



(1.5) 
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where exp{i„4 FV [x + , x_]/7i} is the Feynman- Vernon influence functional. The inluence action i„4 FV [x + ,x_] is the 
sum of a dissipative and a fluctuating part .Af) [x_|_,x_] and yl FV [x + , x_], respectively, whose explicit forms are 



^ v [x + ,x_] = 



2hc 2 



dt / dt' Q(t - t') x + C b (x+ t, x' + t')x' + - x + C b (x+ t, x'_ t')x'_ 



x_C b (x_ i,xi t')x' + - x_C b (x_ i,x'_ i')x'_ 



(1.6) 



and 



A£ v [x+,x_] = 



2hc 2 



dt / dt'Q{t-t') x + A h (x + t,x' + t')x' + -x + A h (x + t,x'_t')x'_ 



x_ A b (x_ t, x' + t')x' + + x A b (x_ t, x'_ i')x'_ 



(1.7) 



where x±, x'_j_ are short for x±(i), x±(t'), and C b (x_ t, x'_ t'), A b (x_ t, x'_ t') are 3x3 commutator and anticom- 
mutator functions of the bath of photons. They are sums of correlation functions over the bath of the oscillators of 
frequency fi k , each contributing with a weight Ck(x)c_k(x') = e lk ( x ^ x ) /2fli i V. Thus we may write 



(x t, x' =J2 c-k(x)c k (x') ( [Xt k (t), Xi(t') 



= —iti 



duj'd 3 k 



{2-kY 
duj'd 3 k 



o- k (Lu')6^e^-*'hmuj'(t-t'), 



(1.8) 



^(xt,x'i')=E c -k( x K( x ')({^ k (t)^ k (t')}) T = ^ / ^^^kW k tr coth 
k 

where cr k (u/) is the spectral density contributed by the oscillator of momentum k: 



TtiJ 
2k B T 



3 i k (x— X ) 



2-7T 

ok(w') = =j?-[<y(w' - fi k ) - + n k )]. 

Zi z k 



coscj'(t-t'), (1.9) 



(1.10) 



At zero temperature, we recognize in (l.S) and (1.9) twice the imaginary and real parts of the Feynman propagator 
of a massless particle for t > t' , which in four- vector notation with k = {lo/c, k) and x — (ct, x) reads 



G(x, x') = - [A(x, x') + C(x, x')] = 



c ik(x-x') m 



(2tt)< 



dujd 3 k ich 



-i\u{t—t')— k(x-x')] 



k 2 + irj J (2tt) 4 uj 2 - nl + it] 



(1.11) 



where r\ is an infinitesimally small number > 0. 

We shall now focus attention upon nonrelativistic systems which are so small that the effects of retardation can be 
neglected. Then we can ignore the x-dependence in (1.8) and (1.9) and find 



C«(xt,x'i') » C**(*,f) = -i£^S^d t S R (t - t'). 
The superscript R indicate that the 5-function is a right-sided one, with the property 

dte(t)S R (t) = 1. 



(1.12) 



(1.13) 



The retarded nature of the <5-function expresses the causality of the dissipation forces, which is crucial for producing 
a probability conserving time evolution of the probability distribution Q . 

Inserting this into (1.6) and integrating by parts, we obtain two contributions. The first is a diverging term 



A„4i oc [x+,x_ 



AM 



dt (x 2 + ~x 2 _)(t), 



where 



AM 



duj'd 3 k er k (a/) x ij 



c 2 J (2tt) 4 u 



i " kk 



6n 2 c 3 



dk. 



(1.14) 



(1.15) 



diverges linearly. This simply renormalizcs the kinetic terms in the path integral (1.5), renormalizing them to 



2 



ft 

By identifying M with M Ten this renormalization may be ignored. 
The second term has the form 



I dt ^n {± 2 + _ k 2_y (L16) 



M 



^ v [x+,x_] =-7yy dt(x + -x_)(*)(*+ +*-)*(*), (1.17) 
with the friction constant 

e 2 2 a , 
birc^M oujm 

where a = e 2 /he ~ 1/137 is the feinstructure constant and u>m = Mc 2 /h the Compton frequency associated with the 
mass M. In contrast to the ordinary friction constant, this has the dimension I/frequency. A dependence of the mass 
renormalization on the p otent ial V (the Lamb shift) will be discussed in Sec. IV. 

The superscript R in (1.17) indicates that the the acceleration (x + + x_)(i) is slightly shifted with respect to the 
velocity factor (x+ — x_)(£) towards an earli er ti me. This accounts for the fact that a physical friction force is always 
retarded, due to the Heaviside function in (1.6). The retardation expresses the causality of the dissipation forces, 
which is crucial for producing a probability conserving time evol ution of the probability distributio n 

We now turn to the anticommutator function. Inserting ( 1.1C| ) and the friction constant 7 from (1.18), it becomes 

e 2 

— A b (xi,x'i') m 2-yk B TK(t,t'), (1.19) 



where 



with a Fourier transform 



K(t,t')=K(t-t')= I — JSVJe"-*^* - ^, (1.20) 

, 2tt 



K{uj i ) ^^L coih ^L i (1.21) 

v ; 2k B T 2k B T' { ' 

whose high-temperature expansion starts out like 

The function K{u>') has the normalization K(0) = 1, giving K(t — t') a unit temporal area: 



dtK(t-t') = l. (1.23) 

Thus K(t — t') may be viewed as a (^-function broadened by quantum fluctuations and relaxation effects. 
With the function K(t,t'), the fluctuation part of the influence functional in (jOI), (1.6), (1-5) becomes 

^V[x + ,x_] = 1^ ^ dt ^ dt' (x+ - x_)(i) K(t, t') (x+ - x_)(i'). (1.24) 

Here we have used the symmetry of the function K (t, t') to remove the Heaviside function Q(t — t') from the integrand, 
extending the range of ^'-integration to the entire interval (t a , tb). We have also introduced the constant 

w = 2Mjk B T, (1.25) 

for brevity. 

At very high temperatures, the time evolution amplitude for the density matrix is given by the path integral 
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£/(x +6 ,x_ b ,t b |x +a ,x_ a ,< a ) = / Px+(i) / Dx_(i) exp 



dt 



M 

T (x^-x 2 _)-(F(x + )-T/(x_)) 



2h 



xexp<-— M7 / dt (x+ - x_)(x + + x_) j / cft(x + -x_) 



2?T 



(1.26) 



where the last term becomes local for high temperatures, since K(t, t') — > 6(t — t'). This is the closed-time path integral 
of a particle in contact with a thermal reservoir. For moderately high temperature, we should include also the first 
correction term in (1.22) which adds to the exponent an additional term. 



24(fc s T) 2 



dt (x + — x_) 



(1.27) 



In the classical limit, the last squeeses the forward and backward paths together. The density matrix ( |l.26j ) becomes 
diagonal. The 7-term, however, remains and describes classical radiation damping. 



II. LANGEVIN EQUATIONS 



For high 7T, the last term in the forward-backward path integral (1.26) makes the size of the fluctuations in the 
difference between the paths y(t) = x+(i) — x_(t) very small. It is then convenient to introduce the average of the 
two paths as x(t) = [x+(t) + x_(t)] /2, and expand 



v(x+|) -^(x-|) ~y W(x)+0(y 3 ) 



keeping only the first term. We further introduce an auxiliary quantity rj(t) by 

f](t) = Mx(t) - M7x(f) + VV(x(i)). 



With this, the exponential function in ( 1.26 ) becomes 

exp 



dty 2 (t) 



(2.1) 



(2.2) 



(2.3) 



where w is the constant (1.25). 

Consider now the diagonal part of the amplitude (2.1 
y& = y<» = 0. It represents a probability distribution 



with Xu 



x_ b = x b and x +a = x_ a = x a , implying that 



P(x & i b |x a t a ) = |(x b ,t b |x a ,< Q )| 2 = J7(x b ,x b ,4|x a ,x a ,t Q ). 



Now the variable y can simply be integrated out in (2.3), and we find the probability distribution 

F[77]cxcxp -i- ^ dtrj 2 (t) 
The expectation value of an arbitrary functional of F[x] can be calculated from the path integral 

(F[x]) v =Af J DxP[»7]F[x], 



(2.4) 



(2.5) 



(2.6) 



where the normalization facto r Af is fixed by the condition ( 1 ) = 1. By a change of integration variables from x(t) 
to r](t), the expectation value (2.6) can be rewritten as a functional integral 



(F[x]) v =M J Vr,P[rj\ F[x], 



(2.7) 



Note that the probability distribution (2.5) is ^-independent. Hence in the approximation ( |2.l| ) we obtain the classical 
Langevin equation. In principle, the integrand contains a factor J _1 [x], where J[x] is the functional Jacobian 



J[x] = Det [6rf(t)/Sx j (t')] = det[(M5 2 - M-/d^ R ) 5 l3 + V i V i F(x(t))]. 



(2.8) 
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It can be sh own that the determinant is unity, due to the retardation of the friction term ||, thus justifying its 
omission in (2.7). 

The path integral (|2.7|) may be interpreted as an expectation value with respect to the solutions of a stochastic 
differential equation (|2.2| ) driven by a Gaussian random noise variable n(t) with a correlation function 



{rf{t)rf (t')) T = 8 l] w 5{t - t'). 



(2.9) 



Since the dissipation carries a third time derivative, the treatment of the initial conditions is nontrivial and will be 
discussed elsewhere. In most physical applications 7 leads to slow decay rates. In this case the simplest procedure to 
solve (2.2) is to write the stochastic equation as 



Mx(i) + W(x(£)) = f)(t) + M 7 x(i), 



(2.10) 



and solve it iteratively, first without the 7-term, inser ting the solution on the right-hand side, and such a procedure 
is equivalent to a perturbative expansion in 7 in Eq. (1.26). 



III. MASTER EQUATION FOR TIME EVOLUTION OF DENSITY MATRIX 

We n ow d erive a Schrodinger-likc differential equation describing the evolution of the density matrix p(x+ a , X-a', t a ) 
in Eq. (1.2). In the standard derivation of such an equation B one first localizes the last term via a quadratic 
completion involving a fluctuating noise variable rj(t). Then one goes over to a canonical formulation of the path 
integral (1.26), by rewriting it as a path integral 



Ur,(x+b,X-b,tb\x+ a ,X-a,ta) = J T>X + (t) J X>X_(*) J J 



x exp 



olt [p+x+ - p x - W r) (p + ,p_,x + ,x_)] 



Then C/^(x + b, x_t,, ib|x +a , x_ a , t a ) satisfies the differential equation 

iWtU^x , y , t\x a , y a , t a ) = H v U v (x ,y , t\x a , y a , t a ). 



(3.1) 



(3.2) 



The same equation is obeyed by the density matrix p(x + , x_; t a ). 

At high temperatures where the action in the path integral ( 1.26| ) is local, we can find directly a Hamiltonian 
without the noise-averaging procedure. However, the standard procedure of finding a canonical formulation is not 
applicable because of the high time derivatives of x(t) in the action of (1.26). They can be transformed into canonical 
momentum variables only by introducing several auxiliary independent variables v = x, b = x, . . . [pH . For small 
dissipation, which we shall consider, it is preferable to proceed in another way by going first to a canonical formulation 
of the quantum system without the effect of electromagnetism, and include the effect of the latter recursively. For 
simplicity, we shall treat only the local limiting form of the last term in (1.26). In this limit, we define a Hamilton- like 
operator as follows: 



P 2 _) + V(x+) - y(x_) + ^2 (x+ - x_)(x + 



. w - 
1 — (x + 



2h 



(3.3) 



Here x, x are abbreviations for the commutators 



x = t[W,x], 



X =-[«,*] 



(3.4) 



A direct differentiation of Eq. (1.26) over time leads to the conclusion that the density matrix p{x+, t a ) satisfies 
the time evolution equation 



ihd t p(x + ,x-;t a ) = Hp{x + ,x-;t a ). 



At moderately high temperatures, we also include a term coming from (1.27) 

w 



U x =i 



24(fc B T) 2 



( X+ -x_r 



(3.5) 



(3.6) 
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For systems with fricition caused by a conventional heat bath of harmonic oscillators as discussed by Caldeira and 
Leggett Hj, the analogous extra term was shown by Diosi ^] to bring the Master e qua tion to the general Lindblad 
form I]] which ensures positivity of the probabilities resulting from the solutions of (3.5). 

It is useful to re-express ( |3.5| ) in the standard quantum-mechanical operator form where the density matrix h as a 
bra-ket representation p(t) — ^2 mn p nm (t)\m) (n\. Let us denote the initial Hamilton operator of the system in (1.1) 
by H — p 2 /2M + V, then Eq. ( |3.5|) with the term (3J3) takes the operator form 



ihdtp = Hp=[H,p\ + — (^XXp + /?XX -X/7X - X/7XJ - y [x, [x,yO|J + 2 ^ kBT y 



x. [x,p]]. 



(3.7) 



The reta rdation of x± in (3.3) leads to the specific operator order in the second term, which ensures that Eqs. (3.3) 
and (3.7) preserve the total probability. 

For a free particle with V = and [H, p] = 0, one has x± = p±/M to all orders in 7, such that the time evolution 
equation (fO) becomes 



ihd t p = [H, p] - 



2M 2 



\P,p]]- 



(3.8) 



In the momentum representation of the density matrix p = J2 PP > Ppp'IpXp'Ij the last term simplifies to — iT = 
—iw(p — p') 2 /2M 2 showing that a free particle does not dissipate energy by radiation, and that the off-diagonal 
matrix elements decay with the rate T. 

In general, Eq. (3.3) is an implicit equation for the Hamilton operator 7i. For small e 2 it can be solved approxiamtely 
in a single iteration step, replacing x by p/M and x = -VF/M in Eq. ([Tt]). 

The validity of this iterative procedure is most easily proven in the time-sliced path integral. The final slice of 
infinitesimal width e reads 



[7(x +b , x_ b , ib|x +OJ x_ a , t b - e) 



dp+(h) f ^P-(4) e i{ p+ (t i; )[x + (t l> )-x + (t 6 -e)]-p_x_-W(t i ,)}_ 



(2tt) 3 



(2tt) 2 



(3.9) 



Consider now a term of the generic form F + (x + )F_(x_) in Tt. When differentiating C/(x + b, x_b, 4| x +ai x -ai 4 — e) 
with respect to the final time tb, the integrand receives a factor —H.(tb). The term F + (x + )F_ (x_) in Ti has the 
explicit form e -1 [F + (x + (tb)) — F + (x + (tb — e))] i^_(x_(tb)). It it can be taken out of the integral, yielding 



e- 1 [F + ( X+ (t b ))U - UF + (x + (t b - e))] F-(x_(t 6 )). 



(3.10) 



In operator language, the amplitude U is equal to U « 1 — ieH/h, such the term F + (x + )i^_(x_) in 7i yields an 
operator 



H,F + ( X+ ) F_(x_) 



(3.11) 



in the differential operator for the time evolution. 

For functions of the second derivative x we have to split off the last two time slices and convert the two intermediate 
integrals over x into operator expressions, which obviously leads to the repeated commutator of TL with x, and so on. 



IV. LINE WIDTH 



Let us apply the master equation (B.7) to atoms, where V(x) is the Coulomb potential, assuming it to be initially 
in an e igenstate \i) of H, with a density matrix p — Since atoms decay rather slowly, we may treat the 7-term 

in (3.7) perturbatively. It leads to a time derivative of the density matrix 



dt(i\p(t)\i) = -J-^W 



|[#,P]PP(0)K) = ^E^IpI/X/IpW = ~ M 7E^ I*. 



fi\ 



(4.1) 



where TtUf i = Ei — Ef and x/i = (/|x|i) are the matrix elements of the dipole operator. A further contribution comes 
the last two terms in (3.7): 



dt(i\p\i) 



w 



<*|p 2 N> + 



\2M 2 {k B Tf 



(# 2 N> 



'fi 



12(fc B T) 2 



(4.2) 



G 



This time dependence is caused by spontaneous emission and induced emission and absorption. To identify these terms 
we rewrite the spectral cont ent o f a single oscillator of frequency in the local approximation (t, t') + (t, t') to 
the correlation functions in (1.8) and (l.E) in the x-independent approximation as 



4/T 



C fe (M') + A b (M') = T ^ 7 (27r)42MfikV ~~~~ 2fcs 



dud k 7r 



( C ° th 2^T + [5{UJ ' ~ Qk) ~ 5{UJ ' + " k)] ^'^'^ 



(4.3) 



[6(u' - n k ) + s(u' + n k )] I e - <w '(*-*'). (4.4) 



as 

n i i \ , , is 4-7T f dud 3 k 7T f . _ . 1 

c b( t, o + A b( t, o = T ^ [8V - no + 

Following Einstein's intuitive interpretation, the first term in curly brackets is due to spontaneous emission, the other 
two terms ac com panie d by the Bose occupatio n fu nction account for induced emission and absorption. In the large-T 
expressions (4.1) and (4.2), the first term in (4.4) leads to part of the sum with ufi > only. This is the famous 
Wigner-Weisskopf formula for the natural line width of atomic levels. 

The rest of the time dependence is due to induced absorption and emission. 



V. LAMB SHIFT 



For atoms, the calculation of the diverging term ( 1.14 ) can be done in a more sensitive approximation. Being 
i nter ested in the time behavior of the pure-state density matrix p — we may calculate the affect of the action 

(1.6) perturbatively as follows: Consider the first term in the local approximation, and integrate the external positions 
in the path integral (11.51) over the initial wave functions, forming 



U, 



ii,ti,;ii : t a 



J dx +b dx_ b J dx +a dx- a (i\x +b ) (i|x_ b )[/(x +b , x_ b , i b |x +a , x_ a , t )(x+&|£) <x_ b |i) 



(5.1) 



To lowest order in 7, the contribution of Q1.6Q can be evaluated as follows 



h 2 c 2 



(5.2) 



Inserting Uu_t a -u,t — e iEi ( tb *> etc., this becomes 



&Uu,t ti *,t a =i^H I " dtdt ' d t dt>C b (t,t'){imt)x'(t')\i} =*iy] / " dtdt' C b (t,t'){i\±(t)\f){f\±'(t')\i). (5.3) 
n c l J ta n c z y J ta 



The integration over t and t' yields 



AU, 



ii,tb;ii,t a ~ ^ 



h 2 c* 



dt 



du'd 3 k 



-ffc 



(-')E 



Ei — Ef — LO 1 — XTj 



(5.4) 



After subtracting the mass renormalization (1.14), which can be rewritten once more in the same form as in (5.4) 
with Ei - E f = 0, the integral yields ln[(A + E t - E f )/\E l - E f \) where A is Bethe's cutoff 0. For A > \E t - E f 
the result (f>A\) imples an energy shift of the atomic level \i): 



A 



(5.5) 



which is the celebrated Lamb shift. Usually, the logarithm is approximated by a weighted average L over energy levels. 
Then contribution of the term (5.4) can be attributed to an extra term in the Hamiltonian (3.7) i?LS = iLMj/2[x, , x]. 
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VI. CONCLUSION 



We have calculated the master equation for the time evolution of the quantum mechanical density matrix describing 
dissipation and decoherence of a point particle interacting with the electromagnetic held. The Hamilton-like evolution 
operator was specified recursively. To lowest order in the electromagnetic coupling strength, we have recovered the 
known Lamb shift and natural line width of atomic levels. In addition, we have calculated the additional broadening 
caused by the coupling of the photons to the thermal bath. 

This research is supported by German HSP III grant. 
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